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1.  Introduction 


An  interim  scientific  report  describing  the  research  activities  supported 
by  the  project  was  submitted  in  July  1982.  A  copy  of  this  report  is  attached 

as  an  appendix  to  this  final  report.  Since  the  interim  report  .was  submitted 

only  a  little  over  three  months  ago  not  much  additional  progress  can  be  re¬ 
ported  at  this  time.  Unfortunately,  both  of  the  graduate  students  whose  work 
was  reported  as  being  in  progress  last  July  have  left  the  project.  Therefore 
the  work  on  Narrow-Band  Adaptive  Array  Processors  (Section  2  of  Appendix  A) 
and  on  the  Frequence-Domain  Implementation  (Section  5  of  Appendix  A)  remains 
incomplete.  The  only  item  on  which  some  progress  has  been  made  has  been  an 
analysis  of  the  performance  degradation  of  adaptive  signal  processors  subject 
to  randomly  varying  changes  in  the  noise  field.  We  include  a  brief  summary 
of  this  work  here.  A  paper  containing  the  results  of  this  phase  of  the 
investigation  is  in  preparation.  [3] 

2.  Adaptive  Array  Processors  in  a  Randomly  Varying  Noise  Field 

The  array  processor  considered  in  the  analysis  is  described  in  Ref.  [2], 

It  consists  of  M  transversal  filters  having  n  weights  each.  Each  of  the  M 

filters  takes  as  its  input  the  signal  from  one  of  the  M  sensors  in  the  array. 
The  outputs  of  all  of  the  filters  are  added  to  produce  the  final  filter  output, 
which  can  be  used  to  either  estimate  target  parameters  or  to  detect  presence 
or  absence  of  a  target  in  a  specified  direction.  In  applications  to  detection, 
the  processor  is  assumed  to  have  information  about  target  direction  and  spatial 
spectrum;  however  it  has  no  information  about  the  noise  to  which  it.  is  supposed 
to  adapt.  In  practice,  this  adaptation  results  in  the  formation  of  nulls  in 


the  direction  of  concentrated  noise  sources. 


The  analysis  is  based  on  the  following  assumption: 

1)  The  received  signal  is  a  white-noise  process;  hence  signals  obtained 
from  adjacent  elements. of  the  delay  lines  used  in  the  transversal  filters 
are  independent. 

2)  The  changing  noise  environment  results  in  a  change  of  the  optimum  filter 
weights.  These  weights  are  assumed  to  be  sample  functions  from  an  independent- 
increment  random  process.  Also,  the  variation  in  the  optimum  weight  are  assumed 
to  be  slow  compared  to  the  adaptation  time  of  the  processor;  this  permits  a 
quasistationary  analysis  of  performance  to  be  used. 

Given  these  assumptions  the  approximate  results  presented  in  our  report 
#  8113  [2]  can  be  shown  to  be  rigorously  correct. 

3.  Papers,  Reports,  Work  in  Progress 

1.  F.B.  Tuteur  and  B.  Heilweil  "Adaptive  Array  Processor  for  Wide-Band  Signals" 
Proceedings  of  the  Workshop  on  Application  of  Adaptive  Systems  Theory,  Yale 
University,  May  1981,  pp.  78-84. 

2.  F.B.  Tuteur,  B.  Heilweil,  and  J.  McMahon  "Adaptive  Array  Processor  with 
Time-Varying  Noise  Patterns",  S&IS  report  No.  8113,  Department  of  Engineering 
and  Applied  Science,  Yale  University,  New  Haven,  Ct.,  September  1981. 

3.  F.B.  Tuteur  "Tracking  Behavior  of  an  Adaptive  Array  Processor  in  a  Randomly 
Varying  Noise  Field".  (in  preparation) 

4.  F.B.  Tuteur  and  J.  McMahon  "Frequency  Analysis  of  an  Adaptive  Array  Processor 
with  Varying  Noise  Field",  (in  preparation) 
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1.  Introduction 

The  last  technical  report  on  research  supported  by  the  project  was 
Technical  Report  No.  S&IS  No.  8113  [1].  This  report  was  submitted  in 
September  of  1981.  In  the  meantime  four  separate  investigations  have  been 
conducted  as  follows: 

1.  The  effect  of  moving  interference  sources  in  narrow-band 
adaptive  array  processors  as  used  in  radar. 

2.  Investigation  of  the  properties  of  adaptive  lattice  filters  with 
a  view  toward  eventual  use  in  adaptive  processors. 

3.  Elimination  of  certain  approximations  in  previous  analyses 

4.  Frequency-Domain  Implementation  of  the  Adaptive  Filter 

A  brief  description  of  these  four  tasks  is  presented  below.  Further  details 
will  be  contained  in  our  final  report,  and  in  several  publications  that  are 
planned . 

2 .  Narrow-Rnrd  Adaptive  Array  Processors 

Task  No.  1  was  an  investigation  of  the  operation  of  the  Howells-Applehaum 
adaptive  algorithm  [2]  when  the  interference  was  moving.  A  preliminary 
conclusion  of  this  investigation  was  that  the  algorithm  as  normally 
implemented  is  quite  insensitive  to  movement  of  the  interference  source,  and 
that  noticable  performance  degradation  would  require  an  interference  source 


moving  past  the  receiver  at  speeds  in  excess  of  March  1  and  at  distances  of  a 
few  feet.  It  appears,  therefore  that  interference  motion  need  not  be 
considered  in  the  design  of  this  type  of  adaptive  loop.  We  hope  to  include  a 
complete  analysis  of  this  problem  in  our  final  report. 

3.  Lattice  Structures 

Lattice  adaptive  filters  have  recently  received  much  attention  in 
applications  such  as  channel  equalization  [3]  and  speech  processing  [4,5]. 

The  lattice  structure  arises  naturally  in  the  theory  of  prediction-error 
filters  [6].  Adaptive  lattice  structures  have  been  proposed  by  Satorius  and 
Alexander  [3],  Griffiths  [7]  and  Hotchkiss  [8],  Their  advantage  over  the 
conventional  transversal-filter  structures  of  Widrow  [10]  are  more  rapid 
convergence,  especially  in  applications  such  as  frequency  tracking  where  there 
may  be  a  large  spread  in  eigenvalues  of  the  signal  covariance  matrix.  This 
superiority  has  been  demonstrated  for  particular  cases  by  computer  simulations 
in  [3,8].  There  would  clearly  be  little  reason  to  consider  optimal  forms  of 
transversal  filters  if  lattice  filters  are  uniformly  superior. 

The  update  equations  for  lattice  filters  are  considerably  more 

complicated  than  those  for  transversal  filters  and  only  very  approximate  and 

preliminary  analyses  of  their  behavior  are  therefore  available.  One  such 

analysis  has  recently  been  published  by  Honig  and  Messer schmitt  [9].  It 

* 

assumes  independence  of  the  various  stages  of  the  lattice  filter  and  makes 
other  simplifying  assumptions  that  are  not  obviously  justified.  However,  one 
conclusion  of  this  work  is  that  the  convergence  speed  of  lattice  structures 
is  not  necessarily  greater  than  that  of  transversal  structures.  Thus  the  issue 
of  the  superiority  of  lattice  structures  appears  to  be  less  clear  cut  then 
some  of  the  optimistic  statements  of  earlier  authors  would  lead  one  to 


believe.  In  any  case  no  analysis  of  lattice  filter  behavior  under  time-varying 
input  conditions  appears  to  have  been  attempted,  and  further  work  on  this  problem 
is  Indicated. 


4.  Improvement  of  certain  derivations 

The  main  result  obtained  in  [1],  Eq.  (53),  is  based  on  a  number  of 
simplifying  assumptions  that  were  not  well  justified.  Specifically,  Eq.  (42) 

is  obtained  on  the  basis  of  the  unproven  assumption  that  cross-product  terms 

— 

are  negligible.  Also  it  is  stated  that  covariance  terms  such  as  X^V^  or  X^e 
are  zero. 

It  has  been  shown  that  these  equations  are  all  exact  if  one  makes  the 
single  assumption  of  a  white  input  spectrum.  More  precisely,  what  is  needed 
is  independence  of  signals  in  adjacent  channels  of  the  delay-line  filter  (ie. 
the  elements  of  the  vector  X).  In  practice  this  is  a  somewhat  restrictive 
assumption  since  one  may  wish  to  consider  colored  inputs.  However,  it  is  an 
assumption  that  has  been  used  in  [10].  Also,  computer  simulations  indicate  that 
the  behavior  of  the  LMS  adaptive  loop  is  not  significantly  altered  if  the  elements 
of  X  are  correlated.  In  any  case,  if  this  assumption  is  made  it  is  easily  shown 
that  past  values  of  the  weight  vector  are  independent  from  current  and  future 
values  of  the  input-signal  vector.  From  this  the  vanishing  of  cross-product  terms 
in  Eq.  (42)  of  [1]  is  easily  demonstrated.  A  further  consequence  of  the  more 
accurate  analysis  is  that  the  condition  ptrR  <  1  may  not  be  sufficient  to  insure 
stability  of  the  filter,  but  that  ptrR  <  1/3  is  definitely  sufficient. 

5.  Frequency  Domain  Implementation 

The  frequency-domain  approach  to  the  basic  IMS  adaptive  algorithm  was 
first  suggested  by  Dentino,  McCool  and  Widrow  [11]  and  nn  analysis  of  this  approach 


was  performed  by  Bershad  and  Feintuch  H2J .  In  a  more  recent  paper  by  Reed  and 
Felntuch  [13]  the  time-domain  and  frequency-domain  approaches  are  compared.  In 
general  the  advantage  of  the  frequency-domain  approach  is  that  it  can  be  taylored 
to  input  signals  that  have  narrow-band  components.  As  far  as  our  specific  problem 
is  concerned,  it  was  hoped  that  the  frequency-domain  approach  might  be  more  satis¬ 
factory  in  dealing  with  the  time-varying  input  covariance  matrix  that  results  when 
interference  motion  is  considered. 


.  The  Freauencv-domain  Adaptive  Arrav  Processor 


The  block  diagram  of  the  frequency-domain  adaptive  array  processor 
is  shown  in  figure  1.  As  in  the  time-domain  filter,  we  have  that  each  of  the 
M  hydrophones  is  connected  to  an  (n+1) -element  tapped  delay-line.  Let  the 
vector  of  signals  in  the  £th  delay-line  at  discrete  time  i  be  denoted  by: 

~L1  =  ^xilO  xill - xiln*  *  1=1 » 2 . M*  ^ 

Each  of  these  M  vectors  is  used  as  the  input  to  a  Discrete  Fourier  Transform 

% 

(DFT),  which  can  be  represented  by  the  transformation: 

DFT 

»  ^ere  (2) 

=  ^iio  ^in  ••••  »  1=1»2,...,M.  (3) 

The  subscript  h  is  used  in  2)  and  3)  because  this  vector  refers  to  DFT 
elements  at  a  specific  hydrophone;  the  vector  defined  in  (3)  will  not  be  used 
in  any  subsequent  analysis.  The  elements  *n  fche  frequency-domain 

processor  are  put  into  vectors  corresponding  to  each  discrete  frequency.  We 
therefore  define  the  vectv 


-im  =  (Xi.lm  Xi2m  XiMm)  » 
The  complex  weight  vectors  are  given  by: 


^i.m  =  (Wilm  WJ2tn  •”  WiMm)  *  m-0,1, . . . ,n. 


The  frequency-domain  adaptive  array  processor  forms  an  estimate  of  the 
n  target  signals  which  are  contained  in  the  delay-line  of  the  hydrophone  that 
is  physically  in  the  middle  of  the  linear  array  (the  delay-line  also  contains 
the  noise,  of  course).  Let  us  assume  that  M  is  odd.  Then  the  processor  forms 
an  estimate  of  the  target  signal  component  of  the  vector  ^(^+1/2)*  '^1*s  *s 

done  by  performing  an  inverse  DFT  on  the  sequence 


Zim  “  ^Lin  ^im  ’  m-0,1‘  <6) 

The  vector  of  values  which  result  from  this  final  inverse  DFT, 


* i  (zi0  zil  •**  Zln) 


(7) 


are  the  estimates  of  the  target  signals. 

It  should  be  noted  here  that  all  the  complex  vectors  defined  so  far 
would,  in  actual  implementation,  be  represented  by  two  vectors,  one  contain¬ 
ing  the  real  part  of  each  element,  and  one  containing  the  imaginary  part. 

The  frequency  domain  weight-update  equation  is  given  by: 

— (i+l)m  "  -im  +  2p (-im  "  -im)  »  m=0»1 . n*  <8> 


In  equation  (8) ,  the  vector  is  equal  to  the  spectral  value  of  the  target 
signal  at  the  discrete  frequency  m  multiplied  by  a  vector  of  phase  shifts 
based  upon  the  known  bearing  of  the  target.  Assuming  again  that  M  is  odd, 
is  given  by 

Pjm  “  S(wm)  (  exp(j  (M-l)Aira/2)  exp(,1  (M-3)Alm/2)  ... 

...  exp(-j  (M-l)  A  /2)  )T  ,  m=0,l . n.  (9) 

im 

where  S(wm)  is  the  spectral  value  of  the  target  signal  at  frequency  m/Ts> 

and 


A.  =  2ir(dmcos0, /c(n+l)T  )  . 

lm  Is 


(10) 


Henceforth  it  will  be  assumed  that  the  target  is  stationary,  so  the 


index  i  will  be  dropped  from  equations  (9)  and  (10)  and  from  in  equation 

(8). 


The  vector  U.  is  given  by: 

— lm 

U,  -  (1/n+l)  Z.  X,  ,  m=0, 1, . . . ,n.  (11) 

— lm  im— lm 

where  the  division  by  (n+1)  appears  because  the  elements  of  are  compared 
in  the  weight  update  equation  to  a  power  spectrum  value. 


5.2.  Analysis  of  Performance  Degradation 

Based  on  analysis  which  will  be  presented  in  a  forthcoming  report,  the 
performance  degradation  of  the  frequency-domain  adaptive  array  processor  ex¬ 
pressed  as  a  fraction  of  the  mean-square  output  is  given  by: 


AJ 

J 


n 

Z 

m=o 


u  +  *  T  -1  *« 

-"TV  tr [P  P  +  R  P  aR  aP  ]  + 
(n+l)  — nt-m _ tn  —m  — m  —m 


(n+l)Ts 

4u 


p  T(R  1r  r  -1)2P  * 
-m  -m  -nun  —m 


1  - 


-/  ,,  s  tr[R  ] 
(n-fl)  — m 


T  -1  * 

Z  P  R  P 
— m  — m  — m 
m=o 


(12) 


where  R  is  defined  as 
—m 

~  )  *  m=0»l . n  ,  (13) 

the  dot  (.)  over  R  in  the  last  term  in  the  numerator  means  time  derivative, 
and  Tg  equals  the  sampling  period,  plots  of  Eg  (12)  have  been  made  for  various 
ratios  of  interference  to  isotropic  noise  power  (c»/l-c»)  and  number  of  array 
elements  (M) .  Signal,  interference,  and  ambient  noise  spectra  were  assumed 


flat  out  to  a  maximum  frequency  f  ;  ie 

m 

Target  spectral  density  S(f)  =  {  Soy/2 

{  0 

Interference  spectral  density  N(f)  =  {  N  /2 

o 

{  0 

Ambient  noise  spectral  density  =  N/2 


I  f  l<  f 
|f|> 


lfl<  £„ 
lfl» 


The  value  of  the  loop  gain  y  was  chosen  to  equate  the  two  numerator  terms 
In  (12),  since  this  is  the  optimum  adjustment.  Under  these  conditions  it 
was  found  that  the  relative  performance  degradation,  plotted  as  a  function 
of  the  quantity  n  T  $>/y  ,  where  y  is  the  effective  beam  width  was  essentially 
the  same  as  in  Fig.  4  of  [1J;  ie.  the  results  for  the  frequency-response 
approach  were  the  same  as  those  for  the  time-domain  approach.  In  view  of  the 
fact  that  similar  assumptions  concerning  spectral  densities  of  the  input 
signal,  etc.  were  used  in  the  two  plots,  this  correspondence  is  perhaps  not 
surprising;  in  fact  it  serves  mainly  to  demonstrate  that  the  two  approaches 
yield  the  same  results. 

In  practice  the  weight-vector  update  for  frequency-domain  adaptive  filters 
cannot  take  place  at  every  sample  instant.  An  update  requires  completion  of  a  DFT 
operation  which  cannot  occur  until  n  time  samples  of  the  signal  to  be  transformed 
have  been  collected  in  a  buffer.  This  extra  delay  tends  to  degrade  the  performance 
of  the  frequency-domain  filter  relative  to  that  of  the  time  domain  filter.  It  is 
easily  shown  that  the  additional  performance  degradation  for  a  constant  velocity 
interference  is  proportional  to  /n  where  n  is  the  number  of  time  samples  needed  to 


fill  the  DFT  buffer. 
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